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MULTIPLICITIES FOR LHC NUCLEAR COLLISIONS 
Qs ! USING HERA STRUCTURE FUNCTIONS 
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We compute in QCD perturbation theory the transverse energy carried by gluons, quarks 
and antiquarks with pt > po ~ 2 GeV in Pb+Pb collisions at t/s = 5500 AGeV by 
using structure functions compatible with the small-x increase observed at HERA. This 
gives a perturbative estimate for the energy and entropy density of the bulk system at 
times r ~ 0.1 fm. The predicted initial gluon entropy density gives a lower limit of 
about 2200... 3400 for the final charged multiplicity. Sources of further entropy increase 
are discussed. 
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The new ep-data from HERA [1-2] imply an increase in the structure function at 
small x relative to the behaviour implied by earlier data. Distribution function analyses 
using the new data have been carried out [3-4]. These analyses, in particular, constrain 
the gluon distribution function, which is not directly measured. 

The purpose of this letter is to estimate, using these new gluon distribution functions, 
the average multiplicity in a central Pb+Pb collision at Large Hadron Collider (LHC) 
energies, y/s = 5500 AGeV. As the first stage, we compute the transverse energy carried 
by minijets, gluons, quarks and antiquarks, with pt > po ~ 2 GeV, produced in these 
collisions [5] . This calculation is as reliable as calculations in lowest order perturbative 
QCD at a scale 2 GeV are: the main uncertainties are higher order corrections, which 
we account for with a K factor, and shadowing, which we neglect. The second stage is 
more model dependent. We assume that the produced gluons form a system in space- 
time at Ti = 1/po = 0.1 fm and calculate its energy density from e{ji) = (E^(\y\ < 
0.5)) x p /(nR\) and estimate from this the entropy density. Assuming that the further 
expansion of the system is adiabatic, we obtain from this a lower limit for the total 
multiplicity: later dynamical phenomena can only increase the entropy. 

The expected charged multiplicity is a crucial variable for the planning of experiments 
[6] and the theoretical predictions lie roughly in the range 1500-8000. Our result with the 
new gluon distributions [3-4] enchanced at small x is dN ch /dy > 2200 for xg(x, ) ~ a; -0-3 
and > 3400 for xg{x) ~ x -0 - 5 ; for the old distributions [7], used in ref. [5], with 
xg(x,4GeV 2 ) ~ const, the limit is only dNch/dy > 910. Note that these estimates 
assume a K factor of 2 and do not take into account shadowing. A reliable error 
estimate would require the study of higher order corrections [8] and gluon shadowing 
[9] for each structure function set separately. These tasks will not be pursued here. 
However, some insight on the uncertainty is provided by the dependence of our results 
on different structure function sets and on the parameter po- I n any case, the new HERA 
distribution functions have a very interesting impact on the LHC heavy ion experiment. 

Proceeding to the computation, all details and the physical motivation are given in [5] . 
Given the distribution functions the basic quantities calculated are the integrated 
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jet cross section 

0jet(Po) = I dp\dy x dy 2 \ 2 ^ jet , (1) 

i p 2 2 dp T dy x dy 2 

and the first £>t moment 

/" dcf) 1 rfcr; et 

£ T )^jet(po) = / dp T dy 1 dy 2 — --: — - (ei+e 2 )PT, (2) 

7 27r 2 dp T dyidy 2 

where e± = e(j>T,yi) = 1 if the parton falls inside the region of our interest and =0 
otherwise. Especially, we are going to consider the central rapidity region \y\ < 0.5. 
The differential cross section is given by 

' ;!i " Yl X lfi/ P ( X liPT) x 2fj/p( x 2,PT)-jf » ( 3 ) 



dp 2 T d yi dy 2 v-,,,^-,™,^ 

9,9.9 



and the distribution functions used are DO', H and D-' by MRS [3-4]. The behaviour 
of the gluon distributions of these sets at the scale Q = pr = 2 GeV is illustrated in 
Fig. 1. 

The results with po — 2 GeV are given in Table 1, and the po dependence is studied in 
Figs. 2 and 3. The numbers of Table 1 are on the pp level, to convert them to numbers 
appropriate for central AA collisions we also need the value of the overlap function 
Taa(^)- We shall use Tp bPb (0) = 32/mb. The estimate for the gluonic contribution to 
initial Et in one unit of y near y = is then 

(E$(\y\ < 0.5)}=T PhPh (0)a jet (p )(E T } 

= 2700 GeV DOl : xg(x) ~ const, 

= 4540 GeV DO' : xg(x) ~ const, (4) 

= 8640 GeV H: xg{x) ~ x" ' 3 , 

= 15300 GeV D— ' : xg(x) ~ x~ ' 5 . 

The increase in E T associated with the small-x increase is significant. However, this 
Et is the final Et only if no further interactions take place in the system. The largest 
decrease in E T is obtained if the system thermalises. 



To obtain a lower limit for the final multiplicity we shall first assume that the system 
thermalises immediately after formation at Ti = 0.1 fm and expands adiabatically there- 
after. Since the entropy is conserved in the comoving frame, we can estimate the final 
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multiplicity from the entropy at the initial time Tj. Aiming at a lower limit we shall only 
consider the dominant gluonic component, which, as seen from Table 1, is about 80% 
of the total. Then pressure, energy density and entropy density are p = aT A , e = 3aT 4 
and s = 4aT 3 with a = 167r 2 /90 = 1.75. Since the longitudinal size of the comoving 
volume is Az = rAy we can estimate e(r;) 
the charged final multiplicity 



(E£(\y\ < 0.5))po/irl£ 2 A and obtain for 



dN ch 2 1 dS 2 1 ~ , „ 

— — « — « TiRlAaTfn 

dy 3 3.6*/ 3 3.6 A 1 



= §3^ ^7r^aTi<Sr(|y| < 0.5)> 
= 914 DOl, 
= 1350 DO', 
= 2180 H, 
= 3360 D -' . 
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The assumption of fast initial thermalization seems reasonable when new HERA struc- 
ture functions are used. In this case the gluonic subsystem is already initially close to 
chemical equilibrium, in the sense that the gluon number density is close to the thermal 
density obtained by assuming that all the initial energy density, e{jj), is thermalised. 
In other words, the average energy of the gluons is close to e/n w 2.7T as in thermal 
gas of bosons. The average energy per gluon within \y\ < 0.5 is obtained from Table 1 
as crjet(-Er(M < 0.5))/2oj e t- The results for different structure functions are 

£/gluon = 4.1 GeV 2.7T, = 2.0 GeV DOl, 
= 3.6 GeV = 2.2 GeV DO', 

(6) 

= 3.2 GeV = 2.6 GeV H, 

= 2.9 GeV =3.0 GeV D-'. 

For the distributions with xg(x) ~ const the thermalisation of gluons requires a degra- 
dation of their average energy by collisions. With the HERA structure functions the 
secondary interactions are only needed for changing the directions of momenta to make 
the distribution uniform. 

The lower limit of multiplicity was obtained by considering only the gluonic component 
and assuming immediate complete thermalisation among gluons. Then the work done 
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against longitudinal expansion is maximal, much of the Et is shifted towards larger 
rapidities and e decreases ~ r~ 4 / 3 . Maximal multiplicity is obtained if the viscosity is 
strong enough to balance the work done by the pressure [10]. In this case all Et remains 
in the interval \y\ < 0.5 and the multiplicity can be estimated by using i?x/particle ps 
0.5 GeV. In a realistic intermediate situation one has to consider, starting from the 
initial conditions in Table 1, the gluon, the quark and the antiquark subsystems and 
their interactions [10,11] which tend to thermalize the system, first the gluons, then 
possibly the quarks and antiquarks. Kinetic and chemical equilibria are attained to 
various degrees. All this is associated with an increase in entropy and multiplicity, but 
quantitative estimates depend on the details of modeling. 

In summary, we have considered the effects of the increase of gluon distributions in the 
small-x region as indicated by the new HERA data on the predictions of perturbative 
QCD - at the small scale of 2 GeV - on the partonic subsystem produced in heavy ion 
collisions at LHC energies. With due reservation concerning possible shadowing effects 
[9], the HERA results give a solid reason for expecting rather large multiplicities and 
thus increase the likelihood of observing new gluon-quark plasma phenomena at LHC. 
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Cjet 
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q 


Q 


D01 


130 
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18.8 


10.0 




12.1 
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1.00 


DO' 
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183 


31.9 
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23.6 
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2.00 
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Viet(E T ) 
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DOl 


994 


787 


136 


70.5 
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84.5 


8.29 


7.41 


DO' 


1570 


1241 


212 


120 




170 


142 


14.5 


13.6 
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3100 


2430 


425 


241 




322 


270 


26.7 


25.5 


D-' 


5890 


4720 


782 


392 




549 


479 


35.8 


34.5 



Table 1. Values of Oj e t(Po) and a- iet {j>o){E T ) in units of mb and mbGeV calculated for 
Po = 2 GeV and for three different MRS distribution functions with xg(x, 4GeV 2 ) ~ 
const (DO'), ~ X" 3 (H) and ~ x~ - 5 (D-') [3-4]. The results with Duke-Owens set 1 
distributions (DOl) [7] are shown for comparison. An overall factor K = 2 is included. 
Both the total as well as the contribution of gluons, quarks and antiquarks is given. For 
each set the upper numbers are for all y, the lower for \y\ < 0.5. 
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Figure Captions 

Fig. 1. Gluon distribution functions MRS DO', D-, and H [3-4] at the scale Q = 2 GeV. 
The old gluon distribution of Duke-Owens set 1 (DOl) [7] is shown for comparison. 

Fig. 2. The integrated jet cross section aj e t(po) computed from eq. (1) with the sets of 
distribution functions marked on the figure and for ^/i = 200 and 5500 GeV. 

Fig. 3. The first pr moment <Xj e t(po) (Et) from eq. (2). The curves are as in Fig. 2. 
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